4 期 应 用 力 学 学 R 
月 Chinese Journal of Applied Mechanics 


Vol. 40 No. 4 
Aug. 2023 


文章 编号 :1000- 4939 ( 2023 ) 04- 0848- 07 


圆柱 约束 下 细 长 压 杆 屈曲 及 后 届 曲 行为 


薄 娜 , 杜 长 城 , 孙 强 强 
(西南 石油 大 学 机 电工 程 学 院 ,610500 成 都 ) 


摘 要 :研究 了 轴 向 力作 用 下 受 圆 柱 横向 约束 的 弹性 杆 的 属 曲 和 后 属 曲 行为 。 通 过 有 限 元 模拟 ,分 
析 了 细 长 压 杆 发 生 正弦 和 螺旋 属 曲 的 轴 向 临界 力 , 提 出 了 正统 和 螺旋 必 曲 临界 点 判定 方法 , 且 与 文 
献 结 果 比 较 验 证 了 分 析 的 正确 性 。 同 时 考察 了 长 细 比 和 边界 条 件 等 因素 对 临界 力 的 影响 。 结 果 表 
明正 弦 属 曲 临界 力 和 螺旋 尾 曲 临界 力 随 杆 的 长 细 比 减 小 而 增 大 。 杆 在 一 定 长 度 范围 内 , 端 部 约束 
条 件 对 临界 属 曲 载荷 及 模 态 的 影响 不 可 忽略 ; 当 杆 的 长 度 足 够 长 时 ,可 以 忽略 边界 条 件 对 临界 载荷 
的 影响 。 
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Buckling and post buckling behavior of slender 


compression rod under cylindrical restraint 


JIANG Na, DU Changcheng, SUN Qiangqiang 


(School of Mechanical Engineering , Southwest Petroleum University ,610500 Chengdu , China) 


Abstract :The buckling and post-buckling behaviors of an elastic rod subjected to both axial force and 
transverse cylindrical constraints are studied. The axial critical forces of sinusoidal buckling and spiral 
buckling of slender rod are analyzed by finite element simulation , and the discrimination method for critical 
point of sinusoidal and spiral buckling is proposed , which is verified by comparison with literature results. 
The influence of slenderness ratio and boundary conditions on the critical force is also investigated. The 
results show that the critical forces of sinusoidal buckling and spiral buckling increase with the decrease of 
the slenderness ratio of the rod. In a certain range of rod length , the influence of end constraint conditions 
on critical load and buckling mode is not negligible. While the influence of boundary conditions on the 
critical load can be ignored as the length of the rod is long enough. 
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横向 圆柱 约束 下 的 弹性 杆 在 工程 应 用 中 普遍 存 
在 ,如 石油 钻井 工程 和 连续 油管 技术 中 用 千 米 长 的 
油管 插 和 人 井 简 "2 用 于 支架 植 人 手术 进入 动脉 血管 
rig eh T 电子 器 件 制 造 中 限制 在 狭 窗 通 道内 的 
纳米 杆 趾 以 及 植物 根系 在 土壤 中 的 生长 中 等 。 横 
向 约束 会 显著 增加 杆 的 轴 向 受 压 承载 能 力 ' ,同时 
导致 压 杆 发 生 二 次 届 曲 出 现 空间 届 曲 构 型 ,上 且 此 时 
压 杆 的 屈曲 行为 受到 多 种 因素 调控 。 因 此 圆柱 约束 
下 的 压 杆 稳 定性 问题 一 直 受 到 国内 外 学 者 的 持续 
关注 。 

LUBINSKI ^ 首次 提出 井 简 中 细 长 杆 的 届 曲 | 
题 ,并 导出 了 其 临界 届 曲 载荷 的 近似 解 。 
DAWSON'”J 研 究 了 和 斜 并 简 中 受 约束 杆 的 正弦 届 曲 和 
螺旋 屈曲 ,得 出 的 屈曲 临界 载荷 公式 得 到 了 较 广 泛 
的 应 用 ,但 公式 中 系数 的 取 值 仍 存在 争议 。 后 续 研 
究 人 员 分 别 用 微分 方程 法 、 梁 柱 模型 和 能 量 法 等 推 
F TRAUE hR. MILLER 4$" 给 出 了 水 平 
圆柱 内 细 长 杆 受 压 的 实验 和 数值 研究 结果 。 近 几 
年 ,LIU 4&'7^ 通过 理论 分 析 、 数 值 模拟 和 实验 测量 
研究 ,提出 了 一 种 新 的 理论 模型 来 分 析 约 束 杆 的 正 
弦 和 螺旋 届 曲 过 程 。XIAO 等 上 1 研究 了 无 摩擦 圆柱 
约束 条 件 下 轴 向 压缩 细 长 弹性 梁 的 三 维 届 曲 和 后 屈 
曲 过 程 ,通过 理论 分 析 和 数值 模拟 得 到 各 构 型 之 间 
过 渡 点 处 的 临界 轴 向 位 移 。 特 征 值 分 析 不 能 考虑 大 
变形 的 影响 ,屈曲 试验 又 往往 受到 设备 精度 和 试 样 
尺寸 的 影响 ;而 采用 有 限 元 法 进行 大 变形 分 析 可 以 
更 准确 的 预测 临界 屈曲 载荷 ,县 能 更 直观 的 显示 杆 
的 力学 行为 , WAE 提出 了 一 种 新 的 边界 
条 件 分 类 方法 ,并 证 明了 第 二 类 情况 下 ,边界 条 件 对 
长 管 柱 的 全 螺旋 届 曲 截面 也 有 影响 。 

本 研究 采用 有 限 元 法 分 析 了 圆柱 约束 下 细 长 压 
杆 在 不 同 边 界 条 件 下 的 届 曲 和 后 届 曲 特性 。 基 于 梁 
单元 对 圆柱 约束 下 的 细 长 压 杆 进行 建 模 , 同时 考虑 
压 杆 变形 时 的 大 位 移 和 接触 非 线 性 ,采用 弧 长 法 进 
行 平 衡 位 移 路 径 搜索 ,考察 了 轴 向 载荷 .长 细 比 和 边 
界 条 件 与 届 曲 和 后 届 曲 的 关系 。 


Z| 


1 力学 模型 


1.1 理论 分 析 


忽略 重力 影响 ,考虑 受 轴 向 力 的 细 长 空心 杆 和 
刚性 圆柱 约束 ,模型 如 图 1 所 示 ( 图 中 虚 圆 环 为 杆 发 


Xem, p. 圆柱 约束 下 细 长 压 杆 屈曲 及 后 屈曲 行为 849 


生 届 曲 后 的 位 置 )。 设 杆 长 度 为 L, 横 截面 面积 为 4， 
抗 弯 刚度 为 B81, 圆柱 约束 的 直径 为 D, 杆 上 端 轴 疝 载 
WA P, 杆 下 端的 轴 向 反 力 为 f.。 杆 发 生 届 曲 后 横 
截面 中 心 的 角 位 移 为 9 , 杆 和 圆柱 约束 中 心 距 离 为 + 
(其 最 大 值 为 杆 与 约束 面 的 间 际 距离 5) , 杆 截面 中 
心 在 XYZ 方向 的 位 移 分 别 为 Uj Ua U3 o 


到 柱 约束 


图 1 杆 和 圆柱 约束 模型 图 
Fig.1 Rod and cylinder constraint model diagram 
用 两 节点 梁 单 元 对 杆 进行 有 限 元 建 模 ,每 个 单 
元 中 每 个 节点 有 6 个 自由 度 , 如 图 2 所 示 。 


一 一 一 一 一 一 一 一 


图 2 有 限 元 分 析 单 元 示意 图 


Fig.2 Schematic diagram of finite element analysis element 


单元 的 总 势能 》 
IT -U-u'F, 
3 il &'odV + [o viv -uF (1) 


其 中 ,U 为 杆 单元 变形 能 。 
T SUE IR ON I 
w = Lu, y Ua; Us; O11» Oi Os; s Uaj Us Uy ,O01 s 07», O31 
(2) 
杆 届 曲 变形 后 ,其 单元 内 的 位 移 量 u 表示 为 
u-N-:u (3) 
其 中 ,N 表示 单元 形 函 数 '“! 。 考 虑 节点 的 大 位 移 ， 
由 横向 位 移 wi \ 双 产生 的 轴 向 压缩 不 可 忽略 ,单元 几 
何方 程 可 表示 为 
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"dg 
Eby Nu 
£p, 
€ — ^ = a =E +E, 
J eA a] 
2 2| dz 
JE 
um 
(4) 


式 中 : sw 是 绕 xz 面 内 弯曲 的 线 应 变 ; sw 是 yz 面 内 


弯曲 的 线 应 变 ; e. 是 轴 向 应 变 ; vy, 是 扭转 前 应 变 ; gp 
ANTH B se, 6; 分 别 为 应 变 的 线性 部 分 和 非 线 性 部 
分 。 单 元 内 应 变 可 表示 为 


= (B, + 3B,)u (5) 


其 中 ,B, A B, 为 几何 函数 矩阵 的 线性 部 分 和 非 线 
性 部 分 。 单 元 的 物理 方程 为 
o =De (6) 
式 中 , D = diag(E,E,E,G) 为 弹性 矩阵 。 将 式 
(2)、(5)、(6) 代 入 式 (1) ,整理 得 到 以 节点 自由 度 
KRIETO K I 。 
杆 与 圆柱 约束 的 接触 判定 及 算法 : 当 杆 的 横向 
位 移 量 yui + wu; =r > 6 时 , 杆 在 xy 横向 与 圆柱 约 
束 面 产生 接触 ,在 有 限 元 分 析 中 忽略 杆 与 接触 面 的 
摩擦, 采用 罚 函 数 法 进行 接触 计算 。 在 势能 泛 函 中 
增加 一 个 惩罚 势能 


I. j? (u)'E,P(u) 


(7) 

其 中 :Es 是 惩罚 因子 ;P RARE, Kfs 
价 于 无 约束 优化 问题 

minl” Cu) = M (u) + I,(u) (8) 

将 单元 节点 位 移 向 量 u^ 中 的 元 素 依次 表示 为 


u(i = 1,,8) , 则 单元 的 市 点 力 为 
Fi = CLEE = 1,.…,8 (9) 
ðU; 


由 此 可 以 组 集结 构 的 节点 力 向 量 从 而 得 到 结构 
的 平衡 方程 组 
F,(u) =0 = 1,2,3,.…,2n 
其 中 :zx 为 结构 届 曲 位 移 问 量 ;n 为 节点 总 数 。 
计算 中 的 一 个 难点 在 于 必须 同时 实现 对 后 屈曲 
平衡 位 移 路 径 的 搜索 并 保证 大 位 移 及 接触 非 线性 计 


(10) 
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算 的 收 傅 性 和 精度 。 本 研究 采用 弧 长 控制 的 Now- 
ton-Raphson 迷 代 法 对 非 线 性 结构 方程 组 进行 计算 。 
其 中 弧 长 法 能 够 实现 对 后 届 曲 位 移 路 径 的 搜索 ,而 
迭代 法 能 够 保证 非 线性 方程 的 求解 精度 。 


1.2 加载 方 式 

条 件 为 两 端 固定 边界 条 件 , 杆 和 约束 的 轴 
线 对 齐 。 加 载 过 程 分 为 3 个 静 力 学 结构 分 析 工 沈 。 
为 了 减 小 载荷 对 杆 的 冲击 力 ,在 载荷 -时 间 历 程 中 采 
用 正 区 波 加 载 方 式 , 如 图 3。 轴 向 力 、 扰 动力 随时 间 


的 函数 表达 式 为 
0,0zx:«1 
FG) 2 4E [sin - 7) +05|,1 <t<3 


F,3ztzxd4 
(11) 


PTE P (12) 
0,2xtx4 
在 杆 中 部 施加 随机 扰动 力 ,扰动 力 的 施加 时 


间 \ 位 置 . 大 小 和 方向 对 结果 几乎 没有 影响 。 


图 3 
Fig.3 Loading mode of F, fo 


F, Vf 加 载 方式 示意 图 


2 屈曲 结果 分 析 


2.1 计算 参数 


设 压 杆 内 径 7 =2 mm, 4$ r, 23 mm, KE L = 
3.5 m, SE PERLE E 2 210 GPa, TRA EC u 20. 3, RH 
梁 单 元 对 压 杆 进行 离散 ,单元 数量 3 500 ,单元 长 度 
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第 4 期 蒋 娜 ， 等 : 圆柱 约束 下 细 长 
1 mm; 采 用 线性 离散 刚性 单元 对 圆柱 约束 进行 建 
模 ,圆柱 约束 直径 D=20 mm ,长度 上 =3.5 m, 


2.2 轴 向 力 A 与 杆 屈曲 关系 


在 轴 向 压力 作用 下 , 杆 会 经 历 3 种 平衡 构 型 : 稳 
定 的 直线 二 维 正弦 形 和 三 维 螺旋 形 。 如 图 4 是 杆 
端 部 位 移 AL SMR FRR. 


500 
400 
300 
E 
~ 200 
100 ; [15.901 
574101 
0 
-03 01 05 09 13 17 2.1 
A L/mm 
KA 端 部 位 移 AL 与 轴 疝 载荷 及 的 关系 


Fig.4 Side displacement AL relations with axial load F, 
由 图 4 可 以 看 出 , 杆 届 曲 的 3 个 阶段 分 别 为 
OA AB 和 BG 段 , 且 判 断 临 界 点 的 方法 分 别 是 位 移 
控制 的 载荷 降 或 载荷 控制 的 位 移 水 平 阶 跃 ,以 下 是 
对 3 个 阶段 的 具体 分 析 。 
第 一 阶段 :图 中 OA 段 杆 发 生 弹性 变形 ,绿色 虚 


3 000 


2 000 
a 
& 
1 000 
SOS 
b NY 
" 
E 
(a) A 点 杆 届 曲 的 投影 图 (b) BAFE h 


压 杆 屈曲 及 后 屈曲 行为 851 
线 斜率 表示 杆 的 纵向 抗 压 刚度 (£4AL) o 

第 二 阶段 :临界 点 A 处 曲线 斜率 出 现 突 跳 , 即 杆 
的 纵向 抗 压 刚度 突然 减 小 , 杆 在 二 维 平 面 内 发 生 届 
曲 失 稳 , 杆 端 部 位 移 AL 增加 速度 变 快 ,直到 杆 与 圆 
柱 约 东 发 生 接触 ,此 时 杆 屈曲 构 型 在 xz 平面 为 正弦 
形 ,如 图 5(a) 。 从 初始 直线 形 到 二 维 正弦 形 的 过 渡 
符合 经 典 屈曲 准则 所) , 即 欧 拉 公 式 所 给 届 曲 临界 压 
力 玉 ,=4mB1P。 在 AB 段 中 , 杆 处 于 正 弱 后 屈曲 
阶段 。 

第 三 阶段 :B 点 轴 向 载荷 F, 保持 不 变 且 位 移出 
现 阶 跃 , 杆 发 生 螺 旋 届 曲 , 且 B. 点 为 第 一 个 螺旋 出 现 
点 。B 点 的 螺旋 届 曲 临界 力 Fo = 115.9 N ,与 文献 
[13,19] 中 二 维 到 三 维 过 渡 时 的 临界 力 吻合 较 好 ， 
文献 [13] 中 下 2 8.18m EI/L? = 120.65N ,本 研究 


结果 比 文献 小 了 3.94% ,表明 了 本 研究 仿真 模型 的 
正确 性 ,误差 的 主要 原因 为 理论 分 析 中 采用 连续 接 
触 假设 和 大 变形 模型 与 小 变形 假设 模型 之 间 的 误 
差 。 随 着 F, 的 增加 ,载荷 位 移 曲线 出 现 类 似 现象 ， 
如 C-G 点 等 , 杆 的 届 曲 螺旋 波形 持续 增加 。 图 5 
(b)、(c) 给 出 了 杆 届 曲 构 型 在 BG 点 时 3 个 面 的 
投影 。 


20, 
O 


i2 
2o, 
的 投影 图 (c) G 点 杆 届 曲 的 投影 图 (x:y:z=1:1:50) 


图 5 杆 届 曲 的 投影 图 


(x:y:z=1:1:50) 


Fig.5 Projection of rod buckling (x: y:z 21:1:50) 


图 6 为 4 个 典型 点 杆 的 正弦 螺旋 屈曲 三 维 构 
型 ,图 中 放大 了 挠 度 来 更 清楚 的 显示 螺旋 形态 ,其 
xz=1:1:50。 可 以 看 出 杆 开始 正弦 失 稳 时 ,与 约 
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束 接触 的 第 一 点 在 约束 中 点 处 ,然后 由 点 接触 变 为 
线 接触 ,接着 发 生 螺 旋 届 曲 , 旦 忽略 重力 的 弹性 杆 屈 
曲 形状 始终 是 对 称 的 。 
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差 值 变 小 ,说 明 螺 旋 届 曲 随 着 杆 长 的 增加 更 易 发 生 。 
在 工程 应 用 中 ,为 了 增加 稳定 性 ,应 选用 长 度 较 小 的 
杆 或 在 长 杆 中 间 增 加 稳定 器 以 减 小 相对 杆 长 。 


2.4 条 件 的 影响 


本 研究 取 12 组 典型 长 度 杆 来 研究 边界 条 件 对 
杆 届 曲 和 后 届 曲 的 影响 ,图 8、9 分 别 表 示 3 种 边界 
条 件 下 的 正弦 屈曲 临界 力 和 螺旋 届 曲 临界 力 。 图 中 
实 线 表示 的 理 华 值 和 虚线 表示 的 仿真 值 吻合 较 好 ， 
误差 均 小 于 8.64% , 


600 端 固定 ， 理 论 值 
= 端 固定 ， nid. 
500 一 一 一 端 固定 一 端 铵 支 ， 理 论 值 
p j] wb, (FUA 
图 6 杆 的 三 维 届 曲 构 型 (x:y:z=1:1:50) F i — PRR, MW 
E Tec PIME 
Fig.6 3D buckling configuration of rod( x: y: z 2 1: 1:50) p 7 
E; 300 
2.3 A€HEE L BT RES Es 
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图 7 368 FH UR ZI R03, IE HN 
与 仿真 值 的 误差 计算 公式 为 
FL — FÌ (e 9 2 3 
Ôi = a Fale) x 100% (13) E ^ Xu 
= 图 8 边界 条 件 与 Fa 的 关系 
8 m Fa — ET x 100% (14) Fig.8 Relationship between FÌ, and boundary conditions 
en2 ^ EM 700 
其 中 : 5. 和 5 分别 表 示 正 弦 届 明和 螺旋 届 曲 的 误 mn —— ame. mem 
38; FF. 分 别 为 正弦 和 螺旋 届 曲 临界 力 理论 值 ， Z 500 uui e crede 
正弦 届 曲 临界 力 为 F, = Am EL ,螺旋 届 曲 临界 RET Ll me, mien 
zs - Yn T JA 
HAF, = 8.187 ENL ; F, (e) FE (e) 分 别 为 正 E 
p 
弦 和 螺旋 临界 力 仿 真 值 。 E 
E -正弦 届 曲 临界 力 ， 仿 真 值 E 
一 一 正弦 届 曲 临界 力 ， 理 论 值 100 


临界 力 ， 仿 真 值 
由 临界 力 ， 
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E 
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图 9 边界 条 件 与 Fo 的 关系 
Fig.9 Relationship between F°, and boundary conditions 


如 图 10 ki he F, 270 N 时 不 同 


400 


Z 300 


X15 KTE PRSE S. BH SS PES EE SC I TED F , 
杆 在 yz 平面 投影 为 三 个 波形 ,而 两 端 固 支 条 件 下 杆 


人 EE I. E yz PEK IAE. AKEE P A 

杆 长 度 L/m 支 最 稳定 ,临界 屈曲 力 最 大 ,两 端 贸 支 情况 下 的 届 曲 

图 7 FPES 临界 力 最 小 。 在 长 细 比 较 小 的 杆 的 应 用 中 ,使 用 两 

Fig.7 Relationship between rod length and buckling 端 固 支 可 以 增加 杆 的 稳定 性 。 图 8.9 中 曲线 旺 下 降 


图 7 中 理论 与 仿真 值 误 差 均 小 于 8.43% ,两 者 趋势 ,当世 接近 零 时 , Fa 和 外, 趋 近 于 无 穷 大 。 当 世 
吻合 较 好 ,由 此 也 验证 了 模型 的 正确 性 。 图 7 显示 接近 无 穷 大 时 ,临界 载荷 六 和 F 趋 近 于 0。 
临界 力 随 杆 长 增加 而 减 小 ,两 组 曲线 随 着 杆 长 度 的 图 11,12 表示 杆 长 和 不 同 边界 条 件 下 届 曲 临界 
增加 越 接 近 , 即 正弦 届 曲 临界 力 与 螺旋 届 曲 临界 力 力 差异 率 的 关系 。 差 异 率 的 计算 公式 为 
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Ai -hs |, 100% (15) 。 届 昌 临界 力 差异 率 10% 的 临界 点 ,A、B 点 所 对 应 的 

1 也 取 较 小 值 , 即 当 杆 长 工大 于 3.1m 时 ,可 以 忽略 两 
其 中 ,4 、 4: 为 同 长 度 不 同 边界 条 件 下 FSI 的 差 值 。  ” 端 边 界 条 件 对 临界 力 的 影响 。 当 杆 长 达到 4.5 m 
图 11 12 中 A 和 B 点 分 别 为 正弦 屈曲 临界 力 和 螺旋 。 时 ,达到 工程 近似 5% 差 异 率 ,如 图 11 点 C 所 示 。 
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图 10 不 同 边界 下 杆 届 曲 模型 (x: y: z=1:1:50) 
Fig. 10 The rod buckling model of different boundaries(x: y: z 2 1: 1: 50) 


"i 一 一 Wit RU PRO E 
Te usen 差异 率 10% 
€ — EGUESN 本 研究 采用 有 限 元 仿真 方法 ,分 析 了 圆柱 约束 
y 下 细 长 压 杆 在 不 同 边界 条 件 下 的 屈曲 和 后 屈曲 特 
W 49 性 ,给 出 了 屈曲 临界 点 判定 方法 。 考 察 了 多 种 因素 
AG.1,1096) 对 届 曲 和 后 届 曲 的 影响 ,具体 结论 如 下 。 
| 5 1) 轴 向 力作 用 下 的 圆柱 约束 杆 会 经 历 3 种 平 
e a 衡 构 形 :直线 构 形 、 二 维 正弦 构 形 和 三 维 螺 旋 构 形 ， 
E ANLE MEE. FEA AA HE E t B) 7 S IRA tb] BA ERT [e o, ER 
图 11 F, 与 差异 率 的 关系 舍 控 制 的 位 移 水 平 阶 跃 。 
Fig.11 Relationship between F*, and difference rate 2) 圆柱 约束 下 的 压 杆 长 细 比 越 小 ,正弦 届 曲 临 
nm 一 一 Wis RR ELSE 界 力 和 螺旋 届 曲 临界 力 越 大 。 压 杆 长 度 相对 较 短 
E i a RUM Pelo 时 , 端 部 边界 约束 条 件 不 同 ,临界 屈曲 载荷 及 模 态 不 
^ NE RR i Tu nid [n] ,边界 条 件 的 影响 不 可 忽略 。 当 压 杆 长 细 比 足够 
大 时 ,边界 条 件 对 临界 载荷 的 影响 小 于 1096 ,可 以 
i 4 忽略 边界 条 件 的 影响 。 
- M 3) RAF KERRI TS HET T e HL AT , 
可 以 较为 准确 地 得 到 圆柱 约束 压 杆 在 轴 向 力作 用 下 
20 B(4.63,1096) 的 届 曲 和 后 届 曲 构 形 ,计算 结果 与 理论 值 及 文献 结 
a A ME NES N 果 相 比 都 具有 很 好 的 一 致 性 。 
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